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Abstract. In the context of the Miiller-Israel-Stewart second-order theory for dissipa- 
tive fluids due to Grad, we analyze the effects of thermal conduction and viscosity in 
heavy ion collisions. We contrast the results to those of the first-order theory due to 
Eckart and to Landau and Lifshitz and to those of perfect (ideal) fluid due to Euler. 
c"^ ■ We study the energy density and entropy density evolution of a pion gas produced in 

the heavy ion collisions. The truncated version of the second-order theory is used to 
find the dissipative quantities. 
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1. Introduction 



Heavy ion collisions such as those at RHIC provide a basic tool to study the properties of 
hot and dense matter produced at high energies. Of particular importance is knowledge 
of the space-time evolution of the matter. This is important when one wants to study the 
transition from hadron to quark and gluon degrees of freedom, a state of matter known as 
QGP (quark-gluon plasma), as predicted by QCD [m. One way to study the dynamics 
of the produced matter is by using fluid dynamics [ g, H]. To probe the non-equilibrium 
properties of the matter we need to include dissipative effects in the fluid dynamic modeling 
of heavy ion reactions. It is known even in non-relativistic dissipative fluids that dissipation 
will have an effect on the observables [ 

The standard theory of dissipative fluid dynamics developed by Eckart [||] and by Lan- 
dau and Lifshitz [ ^] exhibits certain undesirable effects. The resulting transport equations 
of dissipative fluxes lead to parabolic equations for heat conduction and shear diffusion. 
Thus it predicts an infinite speed of propagation for thermal and viscous signals. Many 
applications of dissipative fluid dynamics in relativistic nuclear collisions have used the 
standard theory [ 0]. Because of such undesirable effects, it is then necessary to apply an- 
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other thermodynamic theory of irreversible processes that does not present this anomalous 
behavior. Causal theory of dissipative fluids developed by Miiller [ [l4]] and by Israel and 
Stewart [ [l5j] due to Grad [ [l3|] were developed to remedy some of these undesirable fea- 
tures. The resulting equations of causal theory are hyperbolic in structure and lead to causal 
propagation. Thus causal theory seems to be a good candidate to use instead of the stan- 
dard theory. In the causal-type theories the space of thermodynamic variables is extended 
to include dissipative flows. These dissipative flows (heat flows and viscous pressures) are 
considered as independent variables. The entropy four-flow then depends not only on the 
primary variables (number density, energy density, and pressure) but on these dissipative 
flows as well and its production is semi-positive definite. At equilibrium, the entropy is 
maximum. 

2. Basic Features of Relativistic Dissipative Fluid Dynamics 

The basic formulation of relativistic hydrodynamics can be found in literature [ [| |[ 
[Io[ |TH U|]. We will use the natural units: li — c = k& = 1. The metric tensor we use is 
= diag(+, —,—,—). We consider a simple fluid and no electromagnetic fields. This 
fluid is characterized by, 

N%(x) (particle 4-current) , (1) 

T IJV (x) (energy-momentum tensor) , (2) 

S t '(x) (entropy 4-current) , (3) 

where A — l,...,n for the n conserved net charge currents, such as electric charge, baryon 
number, and strangeness. and T pV represent conserved quantities: 

3^ = , (4) 
d,T" v = . (5) 

The above equations are the local conservation of net charge and energy-momentum. They 
are the equations of motion of a relativistic fluid. There are 4 + n equations and 10 + 4n 
independent unknown functions. The second law of thermodynamics requires 

3pS">0 , (6) 

and it forms the basis for the development of the extended irreversible thermodynamics. 

2.1. Tensor Decomposition 

We now perform a tensor decomposition of N%, T^ v , and with respect to an arbitrary, 
time-like, normalized 4-vector u 11 , u^u^ = 1 . The projection onto the 3-space orthogonal to 
u M is denoted by 



A" v = f v - u"u v = A v " , A" v u v = , A" a Al = A" v , = 3 



(7) 
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The tensor decomposition reads: 

N% = nAif + VZ , (8) 

V" = e«"ii v -pA^ + lW^u^+t^ , (9) 

= su lA + ^ , (10) 



where we have defined 



W = (f + hV» , (11) 
t" v = %^-UA" v , (12) 



and the parenthesis notation is defined by 



= X - (A^ +A V ") . (13) 
In this presentation h is the enthalpy per particle defined by 

(14) 

n 

The dissipative fluxes satisfy the following orthogonality relations: 

M^ = 0, u M <f = 0, u tl W = 0, u M t^ = 0, < = . (15) 

In the local rest frame (LRF) where = (1,0), the quantities appearing in the decom- 
posed tensors have the following meanings: 



n A 




u,N» A 


net density of charge of type A , 


(16) 








net flow of charge of type A , 


(17) 


e 






energy density , 


(18) 


P + Tl 




--A T MV 
3 " v 


pressure , 


(19) 


<f 




u v T vX A» x 


heat flow , 


(20) 






T (psi) 


stress tensor , 


(21) 


s 






entropy density , 


(22) 


<3f 




A V 5 V 


entropy flux . 


(23) 



The angular bracket notation is defined by 



A<" v> = 



I(a^+a^)-1a^a 



(24) 



The space-time derivative decomposes into 

y = u»D + y», «"v„ = , (25) 
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where 

D = uPd^ convective time derivative , (26) 

v /j = A^9 V gradient operator . (27) 

In LRF the two operators attain the given meanings. In this rest frame the projector be- 
comes 

A^ = A^ = diag(0, -1,-1,-1), AC = diag(0, 1,1,1) , (28) 
and the heat four flow has spatial components only (cp = (0, q)). 

3. First Set of Equations: The Conservation Laws 

From now on we will consider one type of charge, namely the net baryon number. With 
the help of the orthogonality properties given in the previous section we obtain the con- 
servation laws. The equation of continuity (baryon conservation), d^N^ = 0, the equation 
of motion (momentum conservation), A^dxT vX = 0, and the equation of energy (energy 
conservation), u^d v T^ = are respectively given by 

Dn = -nv„«' J -V/ + V ) "' J - ( 29 ) 
(e + p + U)Du^ = V A '(/' + n)-Vv^ V +7i'' v L> Mv -[A(;DW'' + 2W^Vv« v) ] ,(30) 
De = -^ + p + Yl)\/ M u fJ + K fjV \/ vUM -\/ M W fJ + 2W fJ D Up . (31) 

There are 5 conservation equations and 14 unknown functions. We need 9 additional equa- 
tions to close the system. So far i/ 1 is arbitrary. It can be chosen to be the particle 4-velocity. 
This is known as the Eckart frame or particle frame. In this frame V = 0. Alternatively 
one can choose it to be the 4-velocity of the energy flow. This is known as the Landau 
and Lifshitz frame or energy frame. In this frame W = 0. Using the fundamental ther- 
modynamic equation of Gibbs the energy balance equation can be written in the following 
convenient form: 

Td lt S f, = o l „nr N -m-d fl (f + (fa ft , (32) 



where 






a a 


= u^dpiia 


(4-acceleration of the fluid) 




= A^Api(3 v M jU -3 jU M v ) 


(vorticity tensor) , 




= A^Aj^dvi^+a^v) 


(expansion tensor) , 
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= A a Pe ap = a a « a 


(volume expansion) , 




= 8 a p--A a p0 


(shear tensor) . 
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4. Second Set of Equations: The Transport Equations 

There are two approaches for finding the transport equations in addition to conservation 
laws. The first one is the phenomenological approach and is based on the second law of 
thermodynamics, that is, the principle of non-decreasing entropy. The standard theories of 
dissipative fluid dynamics [ ^j, ^] assume that the entropy four flow is linear in dissipative 
quantities. Hence they are also referred to as first-order theories. The extended theories 
allow the inclusion of terms that are quadratic in the dissipative quantities, II, <y" and 7i MV . 
They are referred to as second-order theories. Using conservation laws u v d M T MV = and 
= we then obtain the expression for the entropy production d^' from which we 
obtain the transport equations by requiring that the entropy production be positive. The 
second approach uses kinetic theory and is based on Boltzmann moment equations. The 
results of transport equations to be presented here are obtained by using a relativistic Grad's 
14 moment approximation [|l3[], However, one can also use the first Chapman-Enskog ap- 
proximation [ [l^] to find the transport equations in the Eckart theory. Both the phenomeno- 
logical and kinetic theory approaches require that the deviations from local thermodynamic 
equilibrium be small, that is, V ,^,71^ and n are small compared to s,p, and n. 

In the standard Eckart theory one obtains the following set of transport equations for 
the bulk viscous pressure, the heat flow and the shear viscous pressure respectively 

n = -CvX > < 34 ) 

<t = XT(^-Du^=-XnT 2 ^'^) , (35) 
7t" v = 2r)v < "« v> , (36) 

where ^(a,e,n), A,(o,e,n), and r|(o,e,n) are the bulk viscosity, thermal conductivity and 
shear viscosity coefficients. These transport coefficients are required to be positive by the 
second law of thermodynamics 



n 2 7g%y 
t,T XT 2 2y\T 



3^=^-^ + ^>0 ■ (37) 



The resulting equations of motion are parabolic, unstable under perturbations and lead to 
an acausal nature of propagation [ [l5j [16]]. This is a paradox since in special relativity 
the speed of light is finite and all maximum speeds should not be greater than this speed. 
This paradox was first addressed by Cattaneo [ p"8| ]. In the above set of equations, if a 
thermodynamic force is suddenly switched off, then the corresponding flux instantaneously 
vanishes, indicating that a signal propagates through the fluid at infinite speed, violating 
relativistic causality. Even in the non-relativistic case, infinite speeds present a problem, 
since physically we expect the signal speed to be limited by the maximum molecular speed. 
To avoid this paradox Cattaneo introduced ad hoc relaxation terms in the phenomenological 
equations. The resulting equations conform with causality and hyperbolicity requirements. 
The only problem was that a theory developed from first principles was needed. It is from 
these arguments that the extended theory of Miiller, Israel and Stewart was developed. 
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The resulting transport equations derived from extended theories differ from those 
of Eckart-type theories: they contain relaxation terms. These relaxation terms make the 
structure of the resulting equations hyperbolic and thus will conform with causality re- 
quirements. Although it may not be reasonable in some situations, we shall assume here 
for simplicity that there are no viscous/heat couplings (i.e. oco = (Xi = in [ |l5]] ). We 
will also assume that there are no couplings of vorticity and acceleration to the heat and 
shear fluxes. In these approximations we are left with a simple but still causal structure 
of transport equations which has the Maxwell-Cattaneo [ |l8[| form of transport equations. 
The Miiller-Israel-Stewart equations reduce, in the Eckart frame, to 

Tn£>n+n ~ -ce , (38) 

x q D<f + <f ~ %r(^f-Du^j , (39) 
XxDn^ + n^ ~ 2 i\&" , (40) 



where 

Tn = CPo, x q =XT$ h T^=2rip 2 . (41) 

Here the Pa(£, «) are the relaxation coefficients of the dissipative fluxes. As before ^, X, r( 
are the transport coefficients. They involve complicated collision integrals, and they also 
depend on the equation of state. The Xa are the relaxation times. They are sometimes 
taken to be the collision time {x co i w l/(«ov), with a being the cross section and v the 
mean particle speed). In general they are different from the collision time. Here Xa will 
be taken to be the time for dissipative fluxes to relax to their equilibrium values. The 
relaxation terms in the extended theories make it possible for one to study the evolution of 
the dissipative quantities. The transport equations are coupled to the evolution equations for 
number density, energy density and momentum. The evolution equations together with the 
transport equations form a quasi-linear, symmetric and hyperbolic system of 14 first order 
partial differential equations. The system of equations is found to fulfill the requirements of 
causality and hyperbolicity [ [l5|, jig ], One also needs to investigate carefully the conditions 
under which the truncated equations are reasonable. 

5. The equation of state and transport coefficients 

In this presentation we study the energy density and entropy density evolution in the 1+1 
Bjorken hydrodynamic limit [ [Tsfa . We therefore consider equations fiy\\) and (f32|). The 
equation of state is that of a massless pion gas. Thus the pressure is given by p = a T 4 
with a = gi,K 2 /90 where gh = 3 is the number of degrees of freedom, the energy density 
and entropy density are given by 8 = 3 a T 4 and s = 4aJ 3 respectively. From the transport 
equations, the bulk viscous equation does not contribute for massless particles, (^ — > 0) [ 
|p. For the 1+1 dimensional Bjorken-type hydrodynamics the heat term in the energy equa- 
tion will not contribute. Thus we need only the shear viscous pressure for this presentation. 
The energy density evolution equation ([H]) becomes 

d -l = Jl±D- + ± (42) 
dx xx 
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where 



4> 
<I> 

d<t> 
dx 



r\/x 



(perfect fluid) , 
(standard theory) , 

(extended theory) , 



(43) 
(44) 

(45) 



where r| = bT' 1 [ Eo[ and (b = n/$f£) where /„ = 93 MeV is the pion decay constant. 
For massless particles P2 = 3/(4p), and this is used in the expression for x K . The energy 
equation can be solved analytically for the perfect fluid and first order (provided r| is con- 
stant) cases. But since we want r| to depend on temperature or time one can then solve the 
equations numerically or first find the temperature evolution as done in [^1[|. In the case of 
the second-order theory we solve the equations numerically. The proper time evolution of 
energy density is given by 



-4/3 



e(x) = e(t ) 
e(t) = {e(T )-4ri/xo} 



-4/3 



(perfect fluid) , 



-4r\/x (first order) 



(46) 
(47) 



Here To represents the instant at which the expansion starts. In Figs. [I] through ^] we 
show the x dependence of energy density e and entropy density s for the three different 
cases: a perfect fluid, a first-order theory of dissipative fluids and a second-order theory of 
dissipative fluids. In Figs. [l]and [|we show the dependence of the time evolution of energy 
density on the initial time To- In Figs. |]and [i] we show the dependence of time evolution 
of entropy density on the initial time To- In both cases we take the initial initial temperature 
to be 200 MeV. 

First we note that the presence of dissipation makes difference. In the earlier times the 
first-order (parabolic) theory predicts a peak in the energy density and entropy while the 
perfect fluid and second-order (both hyperbolic) predict a monotonic decrease. By looking 
at the dependence of £ and s on To, we see that it is in the early times that the difference is 
most pronounced. In the later times there is an indication that the two dissipative theories 
might converge. Therefore knowledge of different time scales is crucial to determine when 
to apply which theory. In heavy ion reactions, where the dynamics is happening in very 
short times, we need to use the hyperbolic theories of relativistic dissipative fluid dynamics. 

The presence of dissipation in heavy ion reactions will have profound effects on the 
space-time evolution of the system. The freeze-out will be delayed. Temperature and 
energy density decrease slower. Enhancement of entropy production will increase the pro- 
duction of final multiplicity since the two can be related. Since the system takes longer to 
cool this will lead to an enhancement in the production of thermal signals (dileptons and 
photons). 
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Fig. 1. The proper time evolution of energy Fig. 2. The proper time evolution of energy 
density with initial conditions: To = 0.7 density with initial conditions: To = 10 
fm/c and T Q = 200 MeV. fm/c and T = 200 MeV. 



6. Conclusions 

In the early stages of collisions, non-equilibrium effects play a dominant role. Thus, non- 
ideal fluid dynamics must be used to accurately describe the evolution of the system. Be- 
cause of the undesirable features of parabolic theories, it is better to use the hyperbolic the- 
ories. Unlike in the first-order theories, where the transport equations are just the algebraic 
relations between the dissipative fluxes and the thermodynamic forces, in the second-order 
theories the transport equations describe the evolution of the dissipative fluxes from an 
arbitrary initial sate to a final steady-state. The presence of the relaxation terms in second- 
order theories makes the structure of the resulting transport equations hyperbolic and thus 
have well-posed initial value problems. The first challenge faced by the second-order theo- 
ries is the increase in the space of thermodynamic variables which brings new coefficients, 
in addition to transport coefficients, in the theory of non-ideal fluid dynamics. However, 
these new coefficients are determined by the equation of state. Like the primary transport 
coefficients which are constrained by the requirement of second law of thermodynamics, 
they are constrained by the requirements of hyperbolicity and hence causality. Thus, in 
principle, one still needs to know about the transport coefficients and the equation of state 
in order to solve the non-equilibrium fluid dynamics problem. The second challenge in- 
volves solving the equations numerically. Finally one would like to compare the results of 
non-equilibrium fluid dynamics to observables. This will require an effort to solve the full 
system of the resulting equations numerically. 

The consequences of non-ideal fluid dynamics, both first-order (if applicable) and 
second-order were demonstrated here using a simplistic situation. We have seen that dis- 
sipative effects will be important in the early stages of the collision dynamics. A more 
careful study of the effect of the non-ideal fluid dynamics on the observables is therefore 
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Fig. 3. The proper time evolution of en- Fig. 4. The proper time evolution of en- 
tropy density with initial conditions: To = tropy density with initial conditions: To = 



0.7 fm/c and T = 200 MeV. 



1 .0 fm/c and 7b = 200 MeV. 



important. Conversely, measurements of the observables related to thermodynamic quan- 
tities would allow us to determine the importance and strength of dissipative processes in 
heavy-ion collisions. 

Here we have used only the truncated version of the transport equations and a simple 
equation of state. A more realistic situation will require careful analysis of both the trans- 
port coefficients and the equation of state which are employed in the full set of the equations 
(including the terms omitted here). It is then that one may have a better understanding of 
when to use either of these theories in the context of relativistic heavy ion collisions. The 
study of non-ideal or non-equilibrium fluid dynamics will be important for constructing 
hydro-molecular dynamic schemes [ fZ2| ]. A resulting hydro-molecular dynamic scheme 
can then be compared to the studies of multi-fluid dynamics [53j]. 
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